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1  Statement  of  the  problem  studied 


This  research  project  focussed  on  the  enhancement  of  current  vehicle  sim¬ 
ulation  capabilities  of  TARDEC.  Problems  addressed  included  development 
of  efficient  numerical  methods  for  solving  equations  of  motion  with  mg  - 
frequency  oscillations,  and  development  of  methods  and  software  for  sensitiv¬ 
ity  analysis  and  optimal  control  of  engineering  systems  including  multibody 
mechanical  systems. 


2  Summary  of  the  most  important  results 


High-frequency  oscillations  in  nonlinear  ODE/DAE  systems  are  a  problem 
in  vehicle  simulation  because  following  the  oscillations  necessitates  the  use 
of  an  extremely  small  timestep.  However,  many  of  the  oscillations  are  not 
important  for  the  overall  numerical  solution.  We  have  shown  that  methods 
based  on  local  linearization  can  fail  because  the  local  eigenstructure  of  t  e 
problems  oscillates  at  the  high  frequency.  Experiments  have  demonstrate 
that  certain  implicit  methods  combined  with  automatic  stepsize  control  can 
damp  out  the  oscillation  safely,  in  regions  where  its  amplitude  is  too  sma 
to  be  important.  The  usual  stepsize  selection  strategies  must  be  modihe 
so  that  they  are  correct  for  the  limiting  high-index  DAE.  There  is  also  some 
theory  to  support  this  technique  for  the  equations  of  motion.  However,  once 
the  stepsize  is  increased,  problems  with  Newton  iteration  convergence  again 
restrict  the  timestep.  These  problems  are  also  due  to  the  rapidly  chang¬ 
ing  local  eigenstructure.  Some  formulations  of  the  equations  of  motion  are 
more  advantageous  than  others  in  terms  of  Newton  iteration  convergence  for 
highlv  oscillatory  systems.  A  coordinate-split  (CS)  method  has  been  devel¬ 
oped' that,  together  with  a  modified  Newton  (CM)  iteration  is  particularly 
effective.  Numerical  results  for  a  number  of  highly  oscillatory  multibo  v 
systems  demonstrate  that  the  new  method  is  particularly  effective  for  highly 
oscillatory  systems  where  the  oscillation  is  small  and  can  be  damped.  Theory 
was  developed  which  explains  the  Newton  convergence  results. 

A  new  family  of  second-order  methods  called  DAE  a-methods  for  solving  the 
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equations  of  motion  for  flexible  mechanism  dynamics  has  been  developed. 
These  methods,  which  extend  the  alpha  methods  for  ODEs  of  structural 
dynamics  to  highly  nonlinear  DAEs,  possess  numerical  dissipation  that  can 
be  controlled  by  the  user.  Convergence  and  stability  analysis  verify  that  the 
new  methods  introduce  no  additional  oscillations  and  preserve  the  stability  of 
the  original  system.  Convergence  of  the  Newton  iteration  is  achieved  via  the 
coordinate-split  modification  to  the  iteration.  The  new  methods  have  been 
shown  to  be  highly  robust  and  efficient  for  a  variety  of  problems  including  a 
bushing  model  and  flexible  slider-crank. 

A  new  software  package,  DASOPT,  was  developed  for  optimal  control  and 
design  optimization  of  large-scale  differential-algebraic  systems.  These  are 
optimization  problems  with  nonlinear  equality  and  inequality  constraints 
which  include  satisfying  the  DAE  system.  There  are  many  applications, 
including  multibody  systems.  In  DASOPT,  the  time  interval  is  divided  into 
subintervals  and. the  problem  is  solved  via  a  multiple-shooting  type  method. 
This  allows  the  use  of  efficient  DAE  software.  The  derivatives  required  by 
the  optimization  are  computed  via  efficient  DAE  sensitivitvity  software  that 
we  developed. 
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